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Abstract 

BPS monopoles which are periodic in one of the spatial directions correspond, via a 
generalized Nahm transform, to solutions of the Hitchin equations on a cylinder. A 
one-parameter family of solutions of these equations, representing a geodesic in the 2- 
monopole moduli space, is constructed numerically. It corresponds to a slow-motion 
dynamical evolution, in which two parallel monopole chains collide and scatter at 
right angles. 
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1 Introduction 



The idea of understanding the slow-motion dynamics of BPS monopoles by studying 
the geodesies on the moduli space of static monopoles has been extensively exploited 
for monopoles in M 3 — see pQ for a review. For example, much is known about the 
4-dimensional moduli space of centred 2-monopoles in M 3 (the Atiyah-Hitchin space), 
about its geodesies, and hence about the mutual scattering of two monopoles. The 
purpose of this note is to look at an example of the corresponding thing for periodic 
monopoles, ie. monopoles on R 2 x S . (This could be visualized as the transverse 
scattering of two parallel chains of BPS monopoles.) In this case [2], the moduli space 
of centred 2-monopoles is again 4-dimensional, an "ALG gravitational instanton of 
type -Do"; but much less is known about its metric; in particular, the geometry is 
believed to have no continuous symmetries (unlike for the Atiyah-Hitchin metric). 
But it does have discrete isometries, and this enables us to identify some geodesic 
submanifolds and geodesies. We shall see that one such geodesic corresponds to the 
head-on collision of two monopoles, resulting (as usual) in 90° scattering. 

The approach is to use a generalized Nahm transform [3j, which says that periodic 
monopoles correspond to certain solutions of the Hitchin equations. The latter are 
what one gets by reducing the (4-dimensional) self-dual Yang-Mills equations down 
to two dimensions. This system has been studied for some time. Long ago, it was 
noted there there are no finite-energy solutions of the SU(2) Hitchin equations on 
M 2 [U E]. Hitchin's comprehensive paper, investigating the equations on compact 
Riemann surfaces, appeared in 1987 [6J. For non-compact gauge groups such as 
SO(2,l), it was pointed out more recently that finite-energy solutions on M? do exist 

In our case [3], we are interested in the SU(2) Hitchin equations on the cylinder 
1 x S 1 , These equations are analysed in section 2 below; section 3 reviews how 
their solutions may be used to construct periodic monopoles, and defines a natural 
totally-geodesic surface S in the moduli space; and section 4 describes the monopoles 
corresponding to points of S, and in particular a geodesic in S which represents the 
head-on collision of two periodic monopoles. Finally, in section 5, we show how a 
particular explicit solution of the Hitchin equations may be obtained as the TV — > oo 
limit of the Nahm data for a finite chain of N monopoles. 

Except for this special case, the Hitchin-equation solutions are not explicit, and 
have to be obtained numerically. The Nahm transform which generates the monopole 
fields is also implemented numerically, to produce the picture of monopole scattering. 
This is analogous to what was done to obtain examples of iV-monopole scattering in 
M 3 — see [1] for examples. 
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2 The SU(2) Hitchin Equations 

The Hitchin equations may be viewed as a 2-dimensional reduction of the 4-dimensional 
self-dual Yang-Mills equations F\2 = -F34, -F13 = -F42, ^14 = ^23- Assume that 
the gauge potential depends only on the two coordinates (x 1 ,^ 2 ), and write 
s = x 1 + ix 2 and = A3 — LA4. Then the SDYM equations reduce to the Hitchin 
equations 

Ds$ = 0, (1) 
F 12 = (2) 

Here -Dj$ := + L4 5 , $], and $* denotes the complex conjugate transpose of <1>. 
Take the gauge group to be SU(2); the fields $ and A s are trace-free 2x2 complex 
matrices, and are smooth on some surface S having local coordinates (s,s). 

We may partially solve the Hitchin equations as follows. First, note that (JT]) 
implies that det($) is holomorphic, so we write det($) = —H(s), and think of the 
holomorphic function H(s) as being fixed. Next, we can use the gauge freedom to 
diagonalize the gauge field F12 = —2iF s g on £ (as long as there is no global topological 
obstruction — but for the case we are interested in, namely E = lxS 1 , there is no 
such obstruction). So we write F\2 = iBa^, where 03 = diag(l, —1), and where B is 
a real- valued function on E. Let us assume that B is not identically zero (the case 
where B = is easy to solve locally). Then the residual gauge freedom is a local 
0(2), consisting of a local U(l) 

<$>^U- l W, U :=ex.p[iu(s,s)a 3 ], (3) 

plus a reflection. The second Hitchin equation ([2]) now implies that $ has the form 



/ 
9 



where fg = H. The residual gauge freedom acts on the functions / and g by 

/ ^ fe~ 2iu , g^ge 2iu (local U(l)), (4) 
/ l— ^ 9i 9 f (reflection). (5) 

The first Hitchin equation ([T]) is then equivalent to the gauge potential A having 
the form A s = aa^ + a<I>, where a(s, s) and a(s, s) are complex- valued functions, and 
where a given by 

2a = g-%g = -f^dgf. (6) 
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Finally, what remains of ([2]) is equivalent to the two equations 

A [RelogG?//)] = 2(1 +4H 2 )(M 2 - |/| 2 ), (7) 

where A = df + d\ = 4 d s dg is the 2-dimensional Laplacian; and 

f- 1 d s (\f\ 2 a)+g- 1 U\g\ 2 a) = 0. (8) 

From (|6j) it follows that, in order for a to be smooth, the functions / and g must 
have the form 

f = p + e^ 2 , g = n-e-W (9) 

where p± are holomorphic functions with fi+fJ,— = H, and ip = ip(s,s) is a smooth 
function. Using the residual gauge freedom (jlj), we can set ip to be real- valued. Then 
([7]) becomes 

A^ = 2(l + 4|a| 2 )(C+e ,/ '-C-e-^), (10) 
where £± = \p-±\ 2 ; and (JSj) can be re- written as 

e-^/ 2 d s (e^V+a) + e^ 2 ^(e^ M-a) = 0. (11) 

To summarize: given H(s), we first choose a holomorphic splitting = p+P- 
(essentially, the only choice is which zeros of H to include in fj, + , and which to include 
in /i_). Then the Hitchin equations ([HE]) are equivalent to the coupled system (fK)l 
[TT]) for the real function ip and the complex function a. Finally, the functions f,g,a 
(and hence $ and As) are determined by (|9j) and 4a = —dgip. Note that the functions 
ip and a are gauge-invariant, except that ip changes sign under the reflection ((SJ). 

3 Nahm Transform 

In this section, we review how certain solutions of the SU(2) Hitchin equations corre- 
spond, via a generalized Nahm transform, to centred periodic monopoles of charge 2. 
Such a periodic monopole solution consists of an SU(2) gauge field Aj = Aj(x,y,z), 
and a Higgs field $ = &(x,y, z) in the adjoint representation, satisfying 

• Dj<& = —\ejkiFkU where is the gauge field obtained from Ay, 

• <5 and Aj are smooth, and periodic in z with period 2ir; 

• locally in some gauge, <3? and Aj satisfy the boundary conditions (i/7r)(log /o)o"3 
0, A x — * 0, A y — * and A z — (IO/tt)^ — > as p — > oo. (Here p and 9 are polar 
coordinates: x + \y = pe 10 .) 
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A generalized Nahm transform [3] relates such monopoles to certain solutions of 
the SU(2) Hitchin equations on the cylinder ExS 11 . Let r£l, and t with period 1, 
denote the coordinates on this cylinder. The solutions are required to satisfy 

• H(s) = 2 cosh(27rs) — K, for some constant K G C; and 

• F\2 — ¥ as r —* ±oo. 

In the first of these conditions, we could introduce a couple of parameters, by writing 
H{s) = 2C cosh(2-7rs) — K, where C is a complex constant. Then \C\ would determine 
the ratio between the monopole size and the spatial period, while arg(C) would 
determine a spatial orientation. In the monopole picture, these two real parameters 
would show up by a slight alteration in the form of the boundary conditions on $ 
and Aj. For simplicity, both of the parameters are omitted here, but it would be 
straightforward to re-introduce them. It is also worth pointing out that restricting 
to SU(2) rather than U(2) Hitchin fields has the effect of centring the corresponding 
monopoles, ie. removing the translation freedom in x, y, z. 

The details of the transform by which one obtains $ and Aj from $ and Aj are 
as follows. Write £ = x + iy. Let ^+ and \P_ be 2 x 2 matrices satisfying the linear 
system 

2dgV+ + 2A S V + - z^ + + C*- - **- = 0, (12) 
29 s ^_ + 2A S ^^ + z^- + - = 0, (13) 

as well as the normalization condition + = I- Here I is the 2x2 

identity matrix, and ( , ) is the L inner product (e,r> = J^JoQ'Tdtdr. Let * 
be the 4x2 matrix obtained by adjoining ^ + and in other words 



Then $ and Aj are given by 

$ = i(*,r*), Aj = {V,djV). (14) 

The moduli space A4 of centred periodic 2-monopoles is 4-dimensional, and it has 
a natural hyperkahler metric [2]. This space is the same (via the Nahm transform) 
as the moduli space A4 of solutions of the Hitchin equations satisfying the various 
conditions listed above. (Strictly speaking, it has not been proved that the natural 
hyperkahler metrics on A4 and A4 are isometric, but this is expected, on general 
grounds, to be true.) Two of the moduli consist of the real and imaginary parts of 
K, and the remaining two determine the function a. 
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In what follows, we shall concentrate on the surface S in Ai = Ai corresponding 
to a = 0. This surface S is a totally-geodesic submanifold of Ai: to see this, we 
may reason as follows. Note, first, that rotation of the monopole field by ir about the 
z-axis (ie. ( i— > — C or 6 ^ 9 + n) preserves the boundary condition on ($, Aj); the 
change in asymptotic behaviour can be compensated by a gauge transformation. So 
this map is an isometry of Ai. From (|12l [T3j) we see that it corresponds to ty+ i— > ^-i- 
and \£_ i — ^ — together with the map 

$ h-> A 5 i-> A s ip \->ip, an -a. (15) 

So (|15|) . which clearly preserves the Hitchin equations and their boundary conditions, 
is also an isometry of the moduli space Ai. It follows that the surface S corresponding 
to a = is a totally-geodesic submanifold of At; or equivalently that the surface 
corresponding to monopoles invariant under rotations by it about the periodic axis, 
is a totally-geodesic submanifold of A4. 

4 Monopoles and Scattering 

For our subfamily S of solutions with a = 0, the Hitchin equations reduce, in effect, 
to 

= 2(£+e^ -£_e"^). (16) 

In this section, we examine solutions of (|16|) . and describe the corresponding periodic 
monopoles; in particular, a geodesic in S representing a head-on scattering process. 

First, there is the matter of defining /U+(s) and which effectively means 

choosing their zeros. The function H(s) has two zeros, and we shall adopt the 
most obvious choice of allocating one of them to /x+(s) and the other to /i_(s): set 
Ijl± = e ns - X±e~ ns , where A± = {K±\jK 2 - 4) /2. In the square root, take \/K 2 - 4 
to be continuous on the complement of the cut K G (—2,2), with Im(\/K 2 — 4) > 
for Im(K) > 0. 

The boundary condition F\2 — > as r —* ±oo is equivalent to a boundary condi- 
tion on ■0, namely £+e^ — £_e - ^ — > as r — ► ±oo; with our choice of n±{s) this gives 
ip —* as r — > +oo and "0 — > log |A_/A+| as r — > — oo. It is easy to see that, with this 
boundary condition (and fixing K) , any small perturbation tp i— > ip + 5ip of a solution 
?/> of (|16[) has to be trivial. For if 5vp is such a perturbation, with <^/> — > as r — > ±oo, 
then L(5ip) = 0, where L is the strictly-negative operator L = A — 2(£_|_e^ + £_e -1 ^); 
so the only solution is Sip = 0. 

Given that ip{r, t) is to be periodic in t, we need to change gauge in order for / 
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and g to be periodic: instead of the expressions ([9]), we set 



g = n-e 



V'/SgiTrtg-ia; 



(17) 



where lo is some real constant. Notice that if A+ and A_ are interchanged, then 
ip i — > —ip, and / and g are interchanged. This is a gauge transformation ([5]), so the 
corresponding solutions of the Hitchin equations are gauge-equivalent. In particular, 
there is no ambiguity (up to gauge-equivalence) on the cut K S (—2,2). 

If K = ±2, then (|16f) admits the solution ip = 0, and we get explicit solutions 
of the Hitchin equations for which F\2 is identically zero. Conversely, it is straight- 
forward to show that these are the only solutions for which F\2 is identically zero. 
These solutions correspond to taking the explicit Hitchin-equation solution for the 
periodic 1-monopole [HI IB], and recycling it by simply doubling the period. They can 
also be obtained as limits of finite monopole chains, as we show in section 5. 

For K ^ ±2, it is not as easy to find explicit solutions of (|16p . For the discus- 
sion that follows, the equation was solved numerically, for various values of K, by 
minimizing the functional 



Here P = P(r, t) is a fixed function, depending on K, introduced simply in order to 
ensure that the integral converges. Note that (|16|) is the Euler-Lagrange equation 
for E[ip], so any local minimum will be a solution of (|16p . The functional (|18p was 
modelled using spectral methods: Chebyshev with n r grid points for r 6 [—L,L], 
and Fourier with nt grid points for t. The field ip approaches its (finite) boundary 
values very rapidly — for example, ip = o(e _27rr ) as r — > oo — and there is little 
loss of accuracy in taking r to lie in a finite range [— L, L], with L of order unity, 
and with ip attaining its boundary values at r = ±L. The numerical minimization 
was effected using a conjugate-gradient method, for various values of the quantities 
L, n r and nt, and the resulting fields ip compared. For any two such fields ip and 
ip, we typically obtain | \ip — ip\ \ < 0.01 \\ip\\, where || • || denotes the supremum norm 
1 | = max( r t ) |^;| . On the basis of this, we believe that for each K, our numerical 
solution ip is within 1% of the actual solution. 

Given such a solution ip, we then implemented the Nahm transform numerically, 
to obtain the monopole fields (4>, Aj). For each point (x, y, z) of a finite 3-dimensional 
grid, the equations (fT2l[T3l) were solved using a relaxation method, and the monopole 
fields were then computed from (|14|). Since \& — > rapidly as r — * ±oo, it is once 
again reasonable to restrict r to a finite range [— L, L\. The partial derivative in (|14p 
is approximated by a simple finite difference on the grid. As a cross-check, we then 



E 



= f -(djiPf + £+e^ + e-e"^ - P drdt 



(18) 
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sampled the Bogomolny equations Dj& + 2 £ jkiFki = at various points on the grid, 
and specifically near the location of the monopoles, again using finite differences for 
the derivatives. For a sufficiently fine grid (both in r, t and in x, y, z) the accuracy is 
better than 1%, in the sense that \D\Q + F23I < 0.01 l-Di^l (and analogously for the 
D2& and equations). As a further check, this time of the boundary condition, 
we computed the quantity 7r|<3?|/ log 2; on a segment of the x-axis; and verified that, 
as required, it approaches unity as x becomes large. (In the K = —2.3 case, for 
example, the numerical result is that 7r|$|/logx = 0.986 at x = 6.5.) These checks 
allow us to be reasonably confident that our numerical procedures give an accurate 
reflection of the true solution. 

The results of these numerical investigations gives the following picture. For each 
K E C, there is exactly one solution t/j of (|16|) : and solutions for distinct K are not 
gauge-equivalent. So the surface S is diffeomorphic to the plane C, on which K is a 
global coordinate. A solution with \K\ 3> 1 corresponds to well-separated monopoles, 
located at the points x + iy = ±y/—K, z = tt (meaning that these are the points 
in M? x S 1 where the Higgs field is zero). Each monopole is roughly spherical in 
shape (see, for example, the K = ±4 pictures in Figure 1). If, at the other extreme, 
K lies on the segment [—2, 2] of the real line, then the monopoles are located on 
the periodic axis x = y = 0, at z = ir/4 and z = 37r/4. They are elongated in the 
x-direction if — 2 < K < (see the K = — 2 picture in Figure 1), and elongated in 
the y-direction if < K < 2. 

To see the transition between the large- 1 if | and small-|if| regimes, one may ex- 
amine the monopoles corresponding to the one-parameter family K E R. In fact, this 
represents a geodesic in S (as we shall see below), and hence also a geodesic in M.. So 
this family describes a slow- motion dynamical evolution of the system [9j [1] , and pro- 
vides yet another example of 90° scattering following a head-on collision — this time 
between two parallel monopole chains, scattering transversely. See Figure 1, which 
plots the surface |$| 2 := tr($ <!>*) = 0.004, for the monopole solutions corresponding 
to various real values of K. The interpretation in terms of parallel monopole chains, 
is that the chains approach each other in the x-direction; the individual monopoles 
coalesce and then separate to form a single chain of 1-monopoles; and finally the 
monopoles in this chain re-coalesce and then separate to form two parallel chains 
which move apart in the y-direction. 

To see that K G R is a geodesic, we show that the map K 1— > K is an isometry. 
First, note that the isometry of M 2 x S l given by x 1— ► x, y 1— > —y, z ^ 2tt — z 
preserves the field equation and boundary condition for periodic monopoles. So this 
map induces an isometry on M. The corresponding map in the Hitchin-equation 
space is $(r,i) i-> <3?(r, 1 - t)* , A r (r,t) \-> A r (r, 1 - t), A t (r,t) i-> -A t (r,l - t), 
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K=-4 



K=-2.3 



K=-2 




Figure 1: two-monopole scattering — snapshots for K = —4, —2.3, —2, 0, 2.3, 4. 

K K, which preserves the Hitchin equations and its associated constraints. The 
corresponding map on ^± is ^±(r, t; x, y, z) = ^(r, 1 — t; x, —y, 2tt — z). So K i— ► K 
is an isometry of S C Ai, and its fixed set Im {K) = is a geodesic. These monopole 
solutions are invariant under the Klein group D2 , consisting of rotations by tt about 
each of the three coordinate axes. The solution for K = has additional symmetry, 
namely rotation by tt/2 about the z-axis. 

In the asymptotic region \K\ 3> 1, the metric on the moduli space is explicit, and 
simple to write down [2]; but in the interior region, it is likely to be rather compli- 
cated. As a starting-point, we have seen that there is a totally-geodesic surface S 
in the moduli space, representing to periodic 2-monopoles which are invariant un- 
der rotation by ir about the periodic axis; and that these correspond, via the Nahm 
transform, to real- valued functions ip(r,t) satisfying (|16j) . together with appropri- 
ate boundary conditions. One geodesic in S represents a head-on collision of two 
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monopole chains. It should be feasible to extend the procedure used in this note, 
to learn more about the metric and its geodesies, and hence about more general 
scattering processes. 

5 Limits of Finite Chains 

In this section, we will examine how the infinite chain of 1-monopoles arises as the 
N — > oo limit of a finite chain of N single monopoles. One knows that this limit 
is delicate, and that care is needed to avoid divergences (cf. pUl El E]). The aim 
here is to describe the limit in the Nahm-transformed picture. The Nahm data for 
the TV-chain are known [T2], and we will see how to take the N — ► oo limit, thereby 
obtaining the Hitchin data corresponding to an infinite chain. In fact, one may regard 
this monopole field as being a chain of m-monopoles, for any integer rn > 1, with 
U(m)-valued Hitchin data; we shall do the two cases m = 1 and m = 2, and the 
latter therefore links up with our earlier material on 2-monopole chains. 

The Nahm data for a monopole of charge N consists of three anti-Hermitian 
matrix- valued functions T a (u), for a = 1, 2, 3 and u S (—1,1), which satisfy the Nahm 
equation dT a /du = \e a \, c \Tb^ T c ], plus appropriate boundary conditions. The Nahm 
data for an iV-chain have the form T a (u) = —if a (u)J a , where J a are the generators of 
the TV-dimensional irreducible representation of SU(2), and satisfy [J a , Jj,] = ie a bcJc- 
Explicitly, we may set J\ = \{J- + J+), J2 = \^{J- — J+) and {J%)ij = where 

(J+)ij = \5 i+1>j s/{N + 1 + 2i)(N + 1 - 2j), 
(J_)y = iJ ij+ iv / (A r + l + 2j)(iV + l-2i). 

The indices i,j run from — (iV — l)/2 to (N — l)/2, and take values in Z or Z + 1/2 
according to whether N is odd or even. The functions f a are specified in terms of 
Jacobi elliptic functions by 

h{u) = Kk'nc(Ku-k), f 2 (u) = Kk' sc(Ku;k), f 3 (u) = K dc{Ku; k), 

where k is the elliptic modulus, k! = \/l — k 2 the complementary modulus, and 
K{k) the complete elliptic integral of the first kind. These Nahm data correspond to 
a finite-length collinear chain of N single monopoles, and the parameter k determines 
the ratio between the monopole size and the distance between adjacent monopoles 
in the chain. 

Now define a new variable r = Ku/(2n). Then T' a = 2nT a /K solve the rescaled 
Nahm equation 

±T' a = \e ahc [TiT' c \. (19) 
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We restrict attention to the cases where N is odd, and consider the limit N — > oo, 
k' — ► 0, with irNk' — > L for some constant Lei, In this limit the matrices T' a 
become infinite-dimensional, and naturally operate on functions u(r,t), periodic in t, 
with Fourier expansion u = Yljez u j( r ) exp(27riji). Note that the range of r becomes 
(— oo, oo) in the limit. Using known limits of elliptic functions, we obtain 



where A 3 (r,t) = iLsinh(27rr) sin(27r£) and Aj(r, t) = \L cosh(27rr) cos(27rt) . The 
limit of the Nahm equation (|19p can be rewritten as the operator equation 



which is equivalent to the Hitchin equations (pQ G|) under the identifications A\ = 0, 
dt + A2 = limTg, and & = A3 — Explicitly, we obtain the known solution of the 
Hitchin equations corresponding to a charge 1 periodic monopole, namely Ai = 0, 
A 2 = 0, $ = Lcosh(2vrs). 

As mentioned above, one may vary this procedure to obtain Hitchin data for an 
infinite chain of m-monopoles, starting from the same Nahm data; here is the m = 2 
version. This time, define r = Ku/ir, so that T' a = irT a /K solve (fTT?|) . Restrict 
attention to the cases where ./V is even, and consider again the limit where N — > 00, 
k' — > 0, with irNk' — > L constant. The resulting infinite matrices act naturally on 
vectors Uj with j G Z + \ , and hence on functions 



(limT{) • u 
(limT^) • u 
(lim T3) • u 



-Ai(r, t)u{r,t) 

A 3 (r,t)u(r,t), 

d t u(r,t), 



[d r , lim Tq] = \ e abc [lim T' h , lim T' c ) 



(20) 




A direct calculation yields 



(lim 7^) • 
(lim ^2) • 



v 



V 



-A A (r,t)u(r,t), 
A 3 (r, t)u{r,t), 
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Hence we obtain Hitchin data 

n a i7r I 1 \ * L / en[t \ 

^1 = 0, A 2 = , <J> = — cosh(7rs) 

2^0 -1 J 2 K J y e ~ mt J 

With L = 4, these are gauge-equivalent to the explicit solution a = 0, K = —2, 
ip = mentioned in section 4. The above discussion justifies the claim that, in this 
case, the 2-monopole chain is just a 1-monopole chain in disguise. 
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